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Abstract We have observed the nearby isolated globule Barnard 68 (B68) in the (J, K) = (1, 1) and (2,2) inversion lines of 
ammonia. The gas kinetic temperature derived from these is T = 10 ± 1.2 K. The observed line-widths are almost thermal: 
AV = 0.181 ± 0.003 kms~' (AVthetm = 0.164 ± 0.010 kms~'), supporting the earlier hypothesis that B68 is in hydrostatic equi- 
librium. The kinetic temperature is an input parameter to the physical cloud model put forward recently, and we discuss the 
impact of the new value in this context. 
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1. Introduction 



Discovered by Barnard ( 1919| ), the isolated, starless globule 
Barnard 68 (B68) received increased attention recently, after 
Alves et al. ( 2001 ) presented a high resolution extinction map 
and a cloud model, suggesting that B68 has the physical struc- 
ture of a so-called Bonnor-Ebert sphere (BES, i.e. a pressure 
bound, isothermal sphere in hydrostatic equilibrium). With its 
column density and thus number density profiles revealed, B68 
became an id eal object to study molecular depletio n in dark 
core interiors ( Bergin et al. 2002; Hotzel et al. 2002 , hereafter 
Paper I), molecular abundances for a number of species (Di 
Francesco et al. 2002) and the gas-to-dust ratio (Paper I). 

The BES cloud model as based on the measurements of 
|Alves et"aL (2001) fixes the normalised profiles of the den- 
sity («/n c ) and the column density (N/N c ) as functions of the 
normalised radius (r/R), while for the absolute values of the 
central density n c , the central column density N c and the radius 
R additional measurements are necessary. The knowledge on 
any two of the parameters n c , N c , D (distance) and T (kinetic 
temperature) settles the others (Paper I). The column density 
is linked to the extinction profile also via the gas-to-dust ratio 
(we use this term as short form for the more precise hydrogen 
column density per unit reddening by dust). Hence, if the BES 
model holds, one can determine the gas-to-dust ratio if D and 
T are known, or one can determine the distance to the cloud 
by measuring T and applying a "standard" gas-to-dust ratio. In 
any case, the kinetic temperature remains the key parameter to 
measure. A reliable temperature measurement in cold, dense 
cores is possible using the 1.3 cm lines of ammonia (Walmsley 
& Ungerechts |1983| ; panby et al.| |1988[ ). However, previous 
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(Max-Planck-Institut fur Radioastronomie) at Effelsberg. 



ammonia measurements (T = 16 K, Bourke et al. 1995) are 



in co ntradiction to other temperature derivations ( |Avery et al 



1987, Paper I). 



Here we present new ammonia measurements, carried out 
with the Effelsberg 100-m telescope. Apart from the tempera- 
ture derivation and an assessment of the inherent uncertainties, 
we compare the line-width to the width expected from purely 
thermal motion, which is a crucial test of the assumption that 
B68 is in hydrostatic equilibrium. We calculate the distance, 
gas-to-dust ratio and fractional ammonia abundance that fol- 
low from the BES model and the new temperature value. 

2. Observations and data calibration 

The (/, K) = (1, 1) and (2,2) inversion lines of NH3 were ob- 
served simultaneously on May, 6th, 2002 with the Effelsberg 
100-m telecope (40" beam at 23.7 GHz). B68 was between 
12.2° and 15.8°(max) elevation. We used the new 1.3 cm 
HEMT (High Electron Mobility Transistor) receiver in the fre- 
quency switching mode. The system temperature including the 
atmosphere was 63-74 K. The spectrometer was a 8192 chan- 
nel autocorrelator split into four quarters, two for each polariza- 
tion channel, centred on the frequencies of the (J, K) = (1, 1) 
and (2,2) transitions, v (l, 1) = 23.694495487 GHz and 
vo(2, 2) = 23.722633335 GHz respectively (line parame- 
from Kukolich 1967). The velocity resolution was 
kms . Pointing was checked on NGC7027 and 



ters 
0.062 



NRAO530 and was found to be better than 6". 

The data were calibrated using NGC7027, for which we 
assumed an observed flux density of F v (23.71 GHz) = 4.99 Jy 



(Ott et al. 1994). This value takes into account the spatial ex- 



tension of the source (true flux is 1.5 % higher) and the epoch 
of the observation (annual flux decrease is 0.5 %). We multi- 
plied the spectra with a gain-elevation correction factor based 
on the measurements of Altenhoff (1983, private communica- 
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tion), which are consistent with more recent observations at in- 
termediate elevations provided by the telescope team. The val- 
ues for the receiver sensitivity and main beam efficiency were 
0.93 K/Jy and 58 % respectively. The uncertainties of the cal- 
ibration are (1) the uncertainty of the absolute flux density of 
the calibration source, (2) the statistical error of our calibra- 
tion measurements (changes in atmospheric conditions, tele- 
scope gain, etc.) and (3) the remaining uncertainty of the gain- 
elevation dependency at very low elevations. They add up to 
20%. Summing, folding and baseline-fitting were done us- 
ing CLASS. Unfortunately, one of the polarization channels 
showed differences between the on- and off-frequency spec- 
tra (possibly due to a sensitivity gradient of the receiver) and 
also a variation of the spectra with time. As we could not re- 
liably correct for these artifacts, we decided to discard that 
channel in the following analysis. The calibrated spectra, ob- 
tained with a total effective integration time of 134 minutes, 
are shown in Figs, [j] and ||. The beam filling factors assumed 
are 77(1, 1) = 77(2,2)= 1. 

3. Results 

The general concept of deriving cloud parameters from the hy- 
perfine spectra is presented by Ho et al. | ( 1979 ). NH3 fundamen- 
tals can b e foun d e.g. in the papers of Kukolich ( 1967), Poy nter 
& Kakar ( |l975| ) and in the review by |Ho & Townes| ( |l983| ). In 
order to derive reliable values for the intrinsic line-width we 
have convolved the model spectrum with the filter frequency 
response function (approximated by a gaussian with full width 
at half maximum (FWHM) of 1 channel-width) when finding 
the best model spectrum in a ^f 2 -fit to the data. Each channel 
was weighted with the inverse square of the rms (root-mean- 
square) noise of the spectrum determined outside the lines. 
For the (1,1) spectrum, the four unknown variables excitation 
temperature T ex , total opacity Tq(1, 1), line-of-sight velocity Vo 
and line-width AV can be determined from the spectrum. The 
weaker hyperfine lines of the (2, 2) spectrum allow only for the 
fitting of two unknowns, namely tq(2, 2) and Vq. We have as- 
sumed T ex (2, 2) = r ex (l, 1) and AV (2,2) = AV(1, 1) (for a dis- 
cussion of these assumptions see Ho et al.| |1979 ), and pro- 
ceeded modelling the (2,2) spectrum analogously to the (1, 1) 
case. The modelled spectra are shown in Figs. [I]and||, together 
with the observed spectra and the residual spectra. The param- 
eters determined in the ^ 2 -fits are given in Table |jj 

The rotational Temperature T[2, characterising the ratio of 
the populations in the (J, K) = (2,2) and (J,K) = (1,1) rota- 
tional states, can be calculated from 
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Figure 1. B68 as observed in the (J, K) = (1, 1) inversion tran- 
sition of NH 3 . The position observed is a(2000) = 17 h 22 m 46*0, 
6 = -23°49'48". The y-axis gives the brightness temperature, 
7\, . The upper x-axis gives the frequency offset from the central 
frequency of the transition, vo(l, 1), corrected for the relative 
velocity of the cloud. The velocity axis corresponds to vo(l, 1) 
and its zero point refers to vo(l, 1) in the V\ sl - = frame. Top: 
Calibrated spectrum after a total of 134 minutes effective in- 
tegration time. Middle (shifted by -3 K for visibility): ;f 2 -fit 
assuming LTE among the 18 hyperfine components and taking 
the velocity resolution of the spectrometer into acccount. Down 
(shifted by -4 K): Residual spectrum after subtracting the fit. 

Table 1. Primary fitting parameters with formal errors as deter- 
mined in the^ 2 -fits shown in Figs, [l] and [| We label to as the 
total opacity of the transition; it is the the sum of the central 
opacities of the individual hyperfine components (all gaussian 
with FWHM = AV) of the inversion transition in question. The 
small discrepancy between the two transitions in the line-of- 
sight velocities is also present in the (otherwise unused) second 
polarization channel. Possible reasons are a mis-positioning of 
the spectrometers (by 0.2 channels), an inaccuracy in the rest 
frequency difference (of 2.8 kHz or 0.01 %) or a combination 
of rotation (of the globule) and different opacities. 





Spectrum 








(J,*) = (1,1) 


(J,K) 


= (2,2) 


Tex (K) 


6.4 ±0.8 






T (J,K) 


7.2 ± 0.4 


0.22 ± 


0.03 


V a (kms- 1 ) 


3.375 ± 0.002 


3.34 ± 


0.02 


AV (kms" 1 ) 


0.181 ±0.003 







(combining Eqs. (2) and (3) of Ho et al. 1979). We find 
Tu = 9.7 ± 0.3 K. The partition function is calculated assum- 
ing that only metastable {J = K) rotational levels are popu- 
lated and T\2 is characteristic for all metastable levels and 
the ground state. The NH3 column density calculated this 
way is N(NH 3 ) = 5.4 x 10 14 cm 4 . According to Walmsley & (AV^rm = V8 ln2(A:r)/(17;7i H ), where k is the Boltzmann con- 



are almost equal, and we get T = 9.9 K. This temperature cor- 
responds to a thermal line-width of AVtherm = 0.164 km s~' 



Ungerechts (1983) and Danby et al. (1988), T\i can be con- stant and mu the atomic hydrogen mass). See Harju et al 



verted to the kinetic temperature, T. At low temperature they (1993) for the equations used to calculate Af(NH 3 ) and T. 
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Figure 2. B68 as observed in the (/, K) - (2, 2) inversion tran- 
sition of NH3. For observed position, integration time and de- 
scription of plotting axes see Fig. |l| Top: Calibrated spectrum. 
The noise is still too high (and the temperature too low) to see 
the weaker hyperfine components. Middle (shifted by -0.6 K 
for visibility): ;f 2 -fit, the fitted paramters are V and t (2, 2). 
The parameters T ex and AV are kept fixed using the values de- 
rived in the (1, 1) fit. Down (shifted by -0.9 K): Residual spec- 
trum after subtracting the fit. 



The formal error of the temperature reflects only the rms 
noise. The calibration error plays no role, because the two 
lines were measured simultaneously and only their intensity 
ratio enters the temperature derivation. The true uncertainty of 
the kinetic temperature must be assessed through a critical re- 



view of the assumptions involved. As discussed by Ho et al 



(1979), the most serious errors are probably the assumptions 



T ex (l, 1) = ?ex(2, 2) and 7/(1,1) = 77(2,2). As peculiar excita- 
tion conditions are unlikely to prevail in B68 (no internal heat- 
ing sources and no large temperature gradients expected), we 
expect these differences to be < 20 %. Thus we have repeated 
our calculations assuming ±20 % differences between the ex- 
citation temperatures and between the beam filling factors of 
the two transitions. The temperature uncertainty introduced by 
these differences turns out tobe+12/-9%. Including the for- 
mal error from the fit, we finally derive T = 9.9+j'? K. 

We also used the Monte Carlo radiative transfer program 



developed by Juvela (1997) in order to test our previous ap- 



proximation of near-homogeneous excitation conditions along 
the line of sight and this way checking our derived temperature 
va lue. The level energies w ere calculated from the analytical fit 
of Poynter & Kakar ( 1975 ), the radiative rates were ca lculated 
according to formulae in Townes & Schawlow ( 1955 ) and the 
collisional coefficients were taken from Danby et al. ( 1988 ). 

As the absolute density of the underlying cloud model 
scales with the cloud's distance, we have run the program 
with several values of D (between 50 and 300 pc), each time 
optimising the fractional abundance, the turbulent line-width 



and the kinetic temperature. The best fits were obtained for 
T = 9.7 K (independent of D). The lowest^ 2 value was found 
for D = 170 pc, but no accurate determination is possible as 
values between 80 and 250 pc are all consistent with our cali- 
bration uncertainty. We also tried non-isothermal cloud mod- 
els (with unchanged density structure) by coupling the ki- 
netic temperature to the dust temperature (Zucconi et al. 2001) 
in the inner region of the globule (defined by n > n cr i t , try- 
ing « cl -i t = 10 4 , 10 4 5 and 10 5 cm -3 ). These temperature distri- 
butions did not improve the fits. 



4. Discussion 

4. 1. Consistency with other temperature values 

In Paper I we concluded that most likely the kinetic temper- 
ature in B68 is around 8 K. This was based on the measured 
C is O(7 = 1-0) excitation temperature of 8 K and our Monte 
Carlo modelling results. Allowing for the possibility of sub- 
thermal excitation in the less dense outer parts, where most of 
the CO emission comes from, the kinetic temperature derived 
here i s in ag reement with our earlier CO measurements. Avery 
et al. ( 1987| ) deduced from their CO and I3 CO observations an 
outward increasing kinetic temperature between 6 and 1 1 K. As 
they used a constant fractional CO abundance, the derived gra- 
dient must be regarded with caution. However, the given range 
there covers the temperature derived here. Temperatures de- 
rived from ammonia in other globules without internal heating 



sources often lie around 10 K (Lemme et al. 1996). Galli et al 



(2002) calculated the gas temperature distributions in molecu- 
lar cloud cores and predicted for B68 a temperature of about 
10 K, increasing only slightly towards the cloud edge. 

There are two other attempts to derive the kinetic temper- 
ature of B 68 by means of dedi cated NH3 inversion line ob- 
servations: Bourke et al. ( 1995 ) derived T = 16 K. using the 
same underlying assumptions as we have used. As discussed 
in Paper I, this derivation is likely to involve some unfortunate 
error either in the calibration process or in the calculation. Very 
recently, Lai et al. (2002) have presented an estimate of the am- 
monia rotational temperature, which corresponds to a kinetic 
temperature of T = 1 1.2 ± 0.9 K, which is consistent with our 
value. 

4.2. Turbulence 

The thermal line-width AV therm = 0.164 ± 0.010 km s" 1 
is only marginally smaller than the actual line-width 
AV = 0.181 + 0.003 km s _1 , which is an essential condi- 
tion for the BES premise of hydrostatic equilibrium. The 
internal support provided by turbulence is 



E, m h — 



2.33 



AV 
AVL. 



^ therm » 



(2) 



where Etherm is the thermal energy and 2.33f7tH is the mean 
molecular weight we assume throughout this paper. From our 
measurements we derive E tur \, /-Etherm = 3 %. The small contri- 
bution of the macroscopic motion is well below the uncertainty 



4 



S. Hotzel et al.: The kinetic temperature of Barnard 68 



of the thermal energy and therefore it can be neglected as a 
physical parameter of the BES model. 

4.3. Fixing the BES model 

As mentioned in Sect. [j], after determination of the kinetic tem- 
perature the BES model of B68 is fixed if only one of the pa- 
rameters n c , D or gas-to-dust ratio is known. Using Eqs. (2) 
and (1 1) of Paper I we get for T — 10 K 



N c = l.VQxlO^cm-^^) 172 

N c = 2.59xl0 22 cm- 2 ( T ^)- 1 
and 

N - 3 05 x 10 22 crrT 2 ( n ^-)/e(h-k) \ 

l\ c - J.ro X 1U Cm ^ 1.23x10* cm-* mag"' J 



(3) 



In the third equation (the only one which does not require the 
BES assumptions to hold) we have assumed N/N(Hz) = 6/5, 
and the term in brackets corresponds to the standard gas-to- 



dust ratio (Bohlin et al. 1978; Cardelli et al. 1989) 



The number density of molecular hydrogen can be es- 
timated from N H 3 data by balancing collisional excitation 
against emission ( Ho & Townes 1983| ). However, to derive a 
value close to the central density, a high spatial resolution is 
required, which cannot be achieved with single-dish observa- 
tions. As B68 has no detectable foreground stars, its distance is 
only known as far as its association wit h the O phiuchus com- 
plex holds. According to de Geus et al. ( 1989 ) this molecular 
cloud complex extends from D = 80 to 170 pc with a central 
value of 125 pc. The gas-to-dust ratio is expected to vary from 
cloud to cloud, but there are a number of measurements sug- 
gesting that the variation does not exceed a factor of 2 (Paper I). 

The relation between the distance and the gas-to-dust ratio 
imposed on B68 by the BES model reads 



n-xs / N(H 2 )/E(H-K) V 
pC U.23xl0 22 cm- 2 mag-V 



(4) 



The globule is thus located on the near side of the Ophiuchus 
complex unless B68 has a smaller than average gas-to-dust 
ratio. The latter would be unexpected given the common un- 
derstanding of dust evo lution in cold and dense environments 
(Kim & Martin 1996). However, our Monte Carlo simula- 



tions preferred a distance in the 100-200 pc range so this 
possibility should not be ruled out. For the standard gas-to- 
dust ratio the mass and the external pressure of the BES are 
M = 0.9 M (D/85pc) and P R = 2.3 x 10 12 Pa (D/85pc)" 2 
respectively. 

4.4. Ammonia abundance 

The beam averaged ammonia column density has been cal- 
culated in Sect. |5J In order to get the fractional abundance 
^(NH 3 ) = W(NH 3 )/AT(H 2 ), we have convolved the BES col- 
umn density profile with a 40" (FWHM) gaussian. Taking 
into account that the observed position is 19" off the po- 
sition of maximum extinction (as reported by |Bergin et al 
2002), the column density towards the peak extinction position 
is JV(NH 3 ) = 7.7 x 10 14 cirT 2 , and with Eqs. | and | we get 



AT(NH 3 ) = 3.0 x 10~ 8 (D/85pc). This value is almost equal to 
the median value of a sample of 22 ammonia clumps in Orion 
(Harju et al. 1993) and very close to the fractional abundance 
in B217SW, a Taurus dense core with similar temperature and 
size as B68 (iHotzel et al.lkOOlh. 



5. Conclusions 

Our observations support the hypothesis that B68 is in a state 
of isothermal hydrostatic equilibrium. Its kinetic temperature is 
10 + 1 .2 K and its turbulent support is negligible. The ammonia 
abundance is close to the values found in other dark cores, but 
from the BES scaling relations we conclude that either the dis- 
tance or the gas-to-dust ratio of B68 is smaller than previously 
expected. 
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